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Abstract 



We study dynamical surface gravity in a general spherically symmetric setting using Painleve- 
Gullstrand (PG) coordinates. Our analysis includes several definitions that have been proposed in 
the past as well as two new definitions adapted to PG coordinates. Various properties are consid- 
ered, including general covariance, value at extremality, locality and static limit. We illustrate with 
specific examples of "dirty" black holes that even for spacetimes possessing a global timelike Killing 
vector, local definitions of surface gravity can differ substantially from "non-local" ones that require 
an asymptotic normalization condition. Finally, we present numerical calculations of dynamical 
surface gravity for black hole formation via spherically symmetric scalar field collapse. Our results 
highlight the differences between the various definitions in a dynamical setting and provide further 
insight into the distinction between local and non-local definitions of surface gravity. 
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I. INTRODUCTION 



An event horizon is a global property of a black hole spacetime that seems to require, 
in classical general relativity at least, the existence of a singularity It is hoped that the 
singularity inside a black hole will ultimately be resolved in a complete quantum theory of 
gravity, but it is not clear what effect singularity resolution would have on the existence of 
event horizons. For example, the global event horizon may be replaced by a compact trapping 
horizon [IJ [2] . Most work on black hole thermodynamics deals with Killing horizons. These 
are idealizations whose existence is belied in all but the most contrived circumstances by 
the presence of Hawking radiation. While such idealizations are very useful, it is ultimately 
necessary to compute corrections to the thermodynamics that arise from the dynamics of 
the horizon. The thermodynamic properties of black holes therefore need to be formulated 
in a more general setting that allows in principle for the growth by accretion and shrinking 
by radiation of the horizon. One of the first attempts to systematically describe such a 
dynamical setting for black hole event horizons was given in the mid-nineties by Hayward [3] 
who introduced the notion of a "trapping horizon." More recently Ashtekar and collaborators 
formulated a rigorous definition of a "dynamical" horizon |3]. 

Surface gravity is a classical construct that plays an important role in black hole thermo- 
dynamics. It is therefore somewhat surprising that a generally accepted definition of surface 
gravity only exists for Killing horizons in static or stationary spacetimes [HI E] • The standard 
definition is based on the existence of a global time translational Killing vector field that 
becomes null on the event horizon. In a dynamical setting the Killing vector ceases to exist 
so the standard (Killing) definition of surface gravity is no longer valid. Various dynamical 
definitions have been proposed over the years to address this problem [6j 18HT3] but to the 
best of our knowledge no consensus exists on which is correct. A systematic analysis of the 
relative merits of some of these definitions was given by Nielsen and Yoon [B] . The present 
paper builds on this analysis by presenting an analytic and numerical study of various dy- 
namical definitions, including two new ones that have not previously been considered. 

In the first part of the paper we examine the various definitions analytically and classify 
them according to several fairly basic criteria. The criteria we consider are: 

1. Is the definition coordinate invariant? 
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2. Does the proposed definition vanish for extremal horizons? 

3. Is the definition local in the sense of not depending on the matter distribution outside 
the horizon? 

4. Does the definition in general reduce to the standard textbook (i.e. Killing) definition 
for all asymptotically flat black hole spacetimes with a global Killing vector? 

We will see that only two of the five definitions we consider answer 2. in the affirmative. 
Moreover, if the answer to 3. is yes then the answer to 4. must be no. 

The second part of the paper presents a numerical study of dynamical surface gravity 
for horizon formation via the collapse of a spherically symmetric massless scalar field. The 
results demonstrate in a precise physical context the differences between the various def- 
initions in the dynamical region, i.e. while the apparent horizon is growing. In addition, 
the simulation provides a graphic illustration of the qualitative and quantitative difference 
between local and non-local definitions after the horizon has settled down to its final value. 

We work primarily in generalized Painleve-Gullstrand (PG) coordinates, which are one 
representative of a class of coordinates that are regular across future horizons. They are 
distinguished by the fact that spatial slices are everywhere flat. The PG metric takes the 
form: 

. 2 _ ,u ^2 , /\ , l2GM(t,r), _ , 



ds 2 = -a(t, r) 2 dt 2 + ydr + y -^—?-a(r, t)dt j + r 2 dtt 2 , (1) 

where M(t,r) is the generalized Misner-Sharp mass function [Hj. It has a physical inter- 
pretation as the quasi-local energy contained within a sphere of areal radius r at time t. 
In asymptotically Schwarzschild spacetimes the Misner-Sharp mass approaches the ADM 
mass at spatial infinity, whereas in radiating spacetimes it reproduces the Bondi mass at 
future null infinity. The lapse function a(t, r) approaches a function independent of r in the 
asymptotically Schwarzschild region and can be set to one there. This will play an important 
role in the interpretation of the numerical results in what follows. 

The paper is organized as follows: in the next Section we present the definitions of 
dynamical surface gravity that we will be studying and evaluate them in PG coordinates 
in a general spherically symmetric setting. Section III shows that in the context of so- 
called "dirty" black holes (those with non-zero stress energy outside of static horizons) one 
must in general choose between criteria 3. and 4. given above. The extremality condition 
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for dynamical horizons is discussed in Section IV while Section V describes our numerical 
methods and presents the results. The last Section contains a summary and conclusions. 

II. DEFINITIONS OF DYNAMICAL SURFACE GRAVITY 

While some of the following definitions are contained in the papers by Nielsen [TSJ, and 
Nielsen and Yoon [5], we also propose two new definitions that are motivated by our focus 
on PG coordinates. 

A. Static Spacetimes: The Killing Definition 

When the mass function and lapse function depend only on the radial coordinate, the 
spacetime is stationary and possesses a timelike Killing vector, t a . The standard definition 
of surface gravity KKuimg f° r a static black hole with a Killing horizon [7j is: 

t a V 'at 13 = t 13 K Killing, (2) 

where the Killing vector is normalized to have unit norm at spatial infinity. Any static, 
spherically symmetric metric can be written in Schwarzschild coordinates as follows: 

ds 2 = -a(r) 2 f(r)dt 2 s + f^dr 2 + rW, (3) 

where t s is the Schwarzschild time. In the static case it is straightforward to transform to 
PG time t in which the metric takes the form Eq. with GM/r = 1 — / and both M(r) and 
cr(r) functions of r only. Note that the transformation from Schwarzschild to PG time leaves 
the lapse unchanged. As well, we assume asymptotic flatness, so that both M(r) — > m and 
cr(r) — > (Too go to constants at spatial infinity. With the suitably normalized global Killing 
vector t a = (1/ctqo, 0, 0, 0) we have 

where the subscript hor means evaluation on the outer apparent horizon and ' denotes differ- 
entiation with respect to the radial coordinate. If one subscribes to the standard definition 
of surface gravity for static black holes, then all proposed definitions should coincide with Q 
in the limit that the spacetime is static and the horizon is Killing. Without loss of generality 
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(Too can be chosen to be unity and we will henceforth assume that this is the case unless 
otherwise stated. If the spacetime exterior to the horizon is empty, then Einstein's equations 
guarantee that a' = so that a H = = 1. More generally, the ratio (Jul^oo depends on 
the matter content between the horizon and infinity. This is illustrated with some examples 
in Section V. It is for this reason that we refer to the standard definition as non-local. It 
should be noted that Hayward and collaborators [10] have used this non-locality to argue 
against the standard definition and motivate a local definition to be described below. 



B. A Non-Covariant Definition 

The first definition we choose to consider is due to Visser [12] and is defined in terms of 
the PG line element in Q in a slightly different format. Visser writes the PG metric in the 
form 

ds 2 = -c(t, rfdt 2 + (dr - v{t, r)dtf + r 2 dQ 2 . (5) 

In our case, this corresponds to making the substitutions c(t,r) = cr(t,r) and v(t,r) = 
— <r(t, r)yJ2GM{t, r)/r. The outer apparent horizon is then given by c = \v\. Defining a 
quantity g(t) by 

g{t)=\{c 2 - v 2 )\ (6) 

he then argues that the dynamical surface gravity is given by 

By construction this coincides in form with the Killing definition KKuung- Note that here and 
in what follows = implies that all quantities to the right are evaluated on the outer apparent 
horizon. The term "outer horizon" refers to the apparent horizon with the greatest areal 
radius at each spatial slice. 



C. Kodama Vector Definition 

Hayward et al. \X0\ proposed a definition for spherically symmetric spacetimes based 
on the Kodama vector [16J. They argued that this is a natural generalization because the 
Kodama vector field is parallel to the timelike Killing vector field if one exists and under 
certain circumstances coincides with it. The Kodama vector also mimics a feature of the 
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Killing vector that is key to the standard definition: it becomes null on the outer apparent 
horizon. The main disadvantage of this definition is that it is difficult to generalize the 
Kodama vector to non-spherically symmetric spacetimes (although see for example [IT]). 
The Kodama vector is given by 

K» = -^=e^d u r , /i = 0, 1 (8) 

where & v is the volume element associated with the horizon normal directions [18] and 7 is 
the determinant of the 2D horizon normal subspace. In PG coordinates the Kodama vector 
is 

0,0,o). (9) 
Hayward et al. [TU] define the dynamical surface gravity to be: 

«h = -^efdvK* = l -g^V,V v r. (10) 

When evaluated on the outer apparent horizon in PG coordinates this yields 

where the overdot denotes differentiation with respect to the time coordinate. This definition 
is clearly local since it depends only on the mass function and its derivatives at the horizon. 
Although the last term depends on the lapse function it does so in the invariant combination 
adt. 



D. Definitions Using Null Normals 

The next three definitions involve the use of the future/past pointing, ingoing and out- 
going null normals. In PG coordinates, these vectors are calculated to be 



2GM 

F = a 1, -\ a + a, 0, , (12) 



n" = ^l,-y^a-<7,0,oj (13) 

where Z M is the (future) outgoing null vector and is the (future) ingoing null vector for 
a > and /3 > 0. The scaling functions a and are functions of t and r that need to 
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be specified in order to get a unique covariant definition of dynamical surface gravity. The 
cross normalization of the null vectors is: 

Fn,, = -2a 2 a/3. (14) 

We can now also calculate the outward and inward going null expansions in PG coordinates. 
They are, respectively: 



2aa 1 



2GM 



Bi = qTVJv = " " S (15) 



2/3a (l + 



2GM 



9 n = q^V,n u = x — * (16) 

where q^ v is the projector onto the two-surface normal to n M and P. As expected for 
appropriate choices of the signs of a, a and (3 the outward expansion vanishes at r = 2GM 
and is positive if GM/r — > 0, whereas the inward expansion is always negative. We now have 
the machinery at our disposal to examine various covariant definitions of surface gravity. 

The first was proposed by Collins [8] in 1992. He derived a generalized first law that 
related the change of area of an apparent horizon to the change in what he interpreted as 
quasi-local energy. On the basis of this Collins identified the coefficient of the variation of 
the area as a generalized surface gravity, with the result: 

Kc ° l = 6~n ~4GM (1 _ 2GM + — . (17) 
The expression on the far right is valid in PG coordinates and is evaluated on an outer horizon 
r = 2GM. Note that the normalization /3 of the inward going normal has cancelled while 
the normalization of l a remains unspecified. This highlights a main issue in finding a useful 
definition of dynamical surface gravity using global vectors: the choice of parameterization 
for the respective null vectors, which translates to choosing the functions a and /3. Collins 
did not, to the best of our knowledge, address this issue and we will argue below that one 
natural choice is a = 1. 

Fodor et al. [61 [TT] proposed a definition using the outgoing null vector / M as follows: 

Kf p = rvj", (is) 

or alternatively 



In order to specify a and j3, they defined the cross-normalization to be l^n^ = —1 and 
fixed the ingoing null vector by requiring that it be everywhere affmely parameterized and 
normalized in an asymptotically static region such that t^n^ = — 1, where t M is the time 
translational Killing vector normalized to unity at spatial infinity. These conditions result in 
a dynamical surface gravity that takes a particularly simple form in Eddington-Finkelstein 
coordinates, but is difficult to evaluate in the fully dynamical setting in PG coordinates. We 



defer a detailed analysis of this definition to future work. We instead retain ( 19 ) with the 



same cross normalization, but choose a = 1 in (12). The resulting expression for the surface 
gravity is: 

Kp ^^m {1 ~ 2GM "> + 7 (20) 

One feature of this choice is that Z M — > i M as r — > r + (where r + is the outermost apparent 
horizon). Thus the outgoing null normal coincides with d/dt on the outer apparent horizon, 
where t is the PG time, which in turn coincides with the Killing vector in the static case. 

An alternative cross-normalization is provided by Nielsen and Visser [U] who choose 
l^n^ = —2 and the scaling function as a = (3 = (cr) -1 . This is their so called "sym- 
metric" choice of normalization and does not require an asymptotically flat spacetime. The 
expression for the surface gravity in PG coordinates is then given by 

k nv ^^ m {1-2GM'). (21) 

Another definition based on null normals was recently given by Abreu and Visser [T3] . 



The starting point is the same as that of Fodor and Nielsen- Visser (i.e. Eq. (18)) but the 
normalization scheme is different. The outgoing null vectors are fixed in terms of a preferred 
time coordinate, namely the Kodama time. The corresponding expression is most simply 
formulated in Schwarzschild-like coordinates so this definition, like that of Fodor et ai, is 
not well adapted to PG coordinates. What is relevant for our purposes is the interesting 
averaging procedure used by Abreu and Visser. They first defined as the surface gravity 



evaluated via (18), with the same cross normalization as Fodor et al, associated with the 
outgoing future pointing and past pointing null normals respectively. They then took 
the average of these to obtain: 

k av = -(«+ + «-). (22) 



If we start instead from the kpg in Eq. (20) and repeat this averaging procedure, we get a 
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new definition of surface gravity (in PG coordinates): 



a 



AGM 



'1 - 2GM') + 



M 
AM' 



(23) 



Remarkably Eq. (23) is precisely the same as Collins definition (17) with the choice a — 1. 
We will henceforth refer to these two collectively as K nu u for simplicity. Note that the 
expression in (23) differs from that of Hayward's only by a factor of a: K nu u = aKjj. This 



is a crucial difference since it ensures that K nu ii reduces to the Killing definition in the 
stationary limit. 



III. "DIRTY" BLACK HOLES 



By Birkhoff 's theorem, any spherically symmetric solution to the vacuum Einstein equa- 
tions can be put into the form Q with M' = and a' = 0. However, in recent years, 
quantum gravity inspired modifications to Einstein's equations have lead to the consider- 
ation of static, spherically symmetric black hole spacetimes for which there is a non-zero 
effective stress-energy distribution everywhere, including the black hole exterior. The result- 
ing lapse function a = a(r) depends non-trivially on the radial coordinate. These solutions 
are what are referred to as "dirty" black holes [191 ED] • It is straightforward to calculate the 
Killing surface gravity for such solutions by first transforming to PG coordinates and using 
Eq. Q. We will do this in the following for a couple of representative examples. 



A. Stringy Charged Black Hole 



The stringy charged black hole [2T| [22] is a solution to a four-dimensional effective low 
energy theory obtained from string theory. The action contains an electromagnetic field 
and a metric both non-minimally coupled to a scalar dilaton field. The solution we wish to 
consider takes the form: 

ds 2 = -(^- 2 ^dt 2 +(^- 2 ^j \^ + R^R-^Pjdn 2 . (24) 

where Q 2 is the magnetic charge and 0o is the asymptotic value of the dilaton. Note that 



the radial coordinate R in (24) is not the areal radius. To put the metric in PG form, we 



first transform to the areal radius: 



^ = R l R _9^1\ (25 ) 



This yields 



ds 2 = - 1 



2m 

W) 



dt 2 



2m V (d_RV dr2 + r2dn2 



R(r) J \dr J 



where R(r) is: 



and 



R{r) = ^m~ + 



r 2 + 



Q 4 e 4(f>0 
Am 2 



1/2 



dR\ 
dr J 



R(r) (R(r) - 

W) - w 



(26) 



(27) 



(28) 



Eq. (26) is in the form <\3h with 



and 



R(r) J \dr 
fdR\ 

There is one physical horizon located at Rh = 2m, which corresponds to: 



(29) 



(30) 



2m I 2m — 



Q 2 e 2< t>° 



m 



(31) 



Note that the horizon area goes to zero in the limit that Q 2 e 2< ^° = 2m 2 , so that the horizon 
is singular in this limit [22\. 

Using the above formulae in the Killing definition (in fact all global definitions) yields: 



" 1 ft \ 1 

KKitting=^f (r) = — . 

On the other hand, the local Hayward and Nielsen- Visser definitions give: 

i i l - Q 2e20 ° 



(32) 



1 - 



2m 2 



(33) 



As anticipated the two classes of definition disagree when Q 2 ^ 0. Moreover, the local 
definitions give a value for surface gravity that diverges in the singular limit Q 2 e 2 ^ = 2m 2 , 
while the Killing definition yields a finite value. 
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B. Quantum Corrected Black Hole 



This is a non-singular single horizon black hole spacetime derived in [23] using a loop 
quantum gravity inspired quantization scheme similar to those applied to black hole in- 
teriors by a variety of authors [24Tl2~7] . The analytically continued spacetime describes a 
Einstein-Rosen type wormhole whose radius contracts to a minimum value determined by 
the quantum scale before re-expanding inside the black hole interior into an asymptoti- 
cally Kantowski-Sachs type spacetime. This scenario therefore realizes earlier proposals for 
"universe creation" inside black holes [28J. 

The metric, for the quantum corrected black hole in Schwarzschild-like coordinates is 

(34) 

where m is the ADM mass of the black hole while k is a new quantum scale, e = +1 before 
the bounce and e = — 1 after the bounce. It was shown in [23] by a straightforward change 
of coordinates that this analytic continuation produces a metric that is completely regular 
at the bounce. 

In order to calculate the surface gravity we note that for this metric: 
while 

a(r) = (l - - j . (36) 

The horizon radius, 

r H = (Am 2 + k 2 ) 1/2 , (37) 

is always greater than the bounce radius k. One can verify that 

df . 2m 
dr 4m 2 + k 2 ' 

where e = +1 on (before the bounce) and 



(38) 



/4m 2 + k 2 \ 1/2 , . 

• (39) 

The Killing definition yields 

- 1 1 

KKiUin9 ~2{4m 2 + k 2 yi 2 ' ( > 
11 



Hayward's and Nielsen- Visser's definition, on the other hand 

Tfl 

KH = KNV= (4mi + kiy (41) 
In this case both expressions have roughly the same qualitative behaviour, but the local 
form approaches the Schwarzschild value more rapidly as k/m — > 0. 



IV. EXTREMALITY CONDITION IN THE DYNAMICAL SETTING 

One condition that is often imposed is that the surface gravity of extremal black holes 
be zero O [10] . A coordinate invariant definition of extremality that applies to dynamical 
horizons can be given as follows |32j : 

r^V^, = 0. (42) 

Geometrically, a spherically symmetric dynamical horizon is extremal when the correspond- 
ing trapping surface is tangent to the inward going null vector. This in turn guarantees that 
the rate of change of the outward expansion along the inward null vector vanishes. Such 
extremal horizons can in principle form in dynamical settings [29] . 
In generalized PG coordinates a straightforward calculation yields: 

^-^(™ + ^), (43) 

so that dynamical extremal horizons satisfy: 

1 - 2GM' + — = 0. (44) 
a 

Clearly only the Hayward definition and K nu n yield a surface gravity that generically vanishes 
for dynamical extremal horizons. 

One can also consider a slicing dependent notion of "extremality" as follows. For a 
spacelike slicing of an evolving black hole spacetime that is regular across future horizons, 
a given spatial slice can intersect the trapping horizon multiple times (See Figure ([2])). The 
instant of formation is defined as the time when the first spatial slice touches the trapping 
horizon. Future spatial slices necessarily contain an inner and an outer horizon. The former 
necessarily moves inward while the latter expands. At the instant of formation the inner 
and outer horizons coincide and in this sense the evolving black hole can be thought of as 
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extremal. The condition for this to occur is that the trapping horizon be tangent to the 



spatial slice In PG coordinates this can easily be worked out using Eq. (15) to yield the 
condition 

1 - 2GM' = 0. (45) 

Note that in the case of multiple intersections, this condition is also satisfied on formation of 
any additional pair of horizons. Clearly, the definitions by Visser and Nielsen- Visser are the 
only ones that will generically vanish on formation in PG coordinates. This will be verified 
in the numerical simulations below. 



V. NUMERICAL METHOD AND RESULTS 

We consider black hole formation via the collapse of a minimally coupled massless scalar 
field ip(r,t). The action is given by 

S = j d A x^-gR{g) - ± J d 4 x\W^\ 2 . (46) 

As derived in detail in [30J, the evolution equations for the scalar field i/j and its canonical 
conjugate 11^ in PG coordinates (JT]) are: 



, ,2GM „ II,. . 

* = °\\l— ^ + tt)> ( 47 ) 



n 



a | r 2 ^+xj——U lp 



(48) 



where we have absorbed a factor of y/8irG into ip in order to make it dimensionless. 
The mass function M is determined by the Hamiltonian constraint: 



1 (% , 2/ //^2^ — l 2GM 



M' = -^ + rWJ+^'n^^— . (49) 

With a slight abuse of terminology we henceforth refer to M' as the "mass density." Con- 
sistency of the evolution equations in this gauge forces the lapse function to satisfy the 
constraint: 

a H a = 0. 50 

V2GMr 

The above equations are evolved numerically with a slightly modified version of the code 
first used in [30J. 
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The spatial grid resolution does not change with time, but does vary across the lattice. 
Near the origin, where high accuracy is needed, we use a fine mesh with a grid resolution of 
10 -5 . From just outside the origin to beyond the outermost apparent horizon the resolution 
is 10 -4 . Finally, from outside the outermost horizon to the end of the lattice we use a grid 
resolution of 10~ 2 which allows us to extend the lattice far enough to contain all the matter 
while reducing the CPU time for each simulation. This coarse resolution is still accurate 
enough to handle the dynamics and keep the code stable. A smoothing function is applied 
where the resolution changes. 

The time steps are determined using an adaptive refinement method. This method scans 
the spatial slice for each time step and refines the time spacing according to the condition 



where dt/dr is the inverse of the local speed of either the ingoing or the outgoing null 
geodesic (whichever yields a smaller At). The time and spatial integrations are done using 
fourth order Runge-Kutta methods. For the run that we illustrate, there are 3000 time steps 
between the formation of the first set of apparent horizons (t = 1.66) and the formation of 
the singularity (t = 1.87). The code becomes unstable just before singularity formation, so 
the last 100 or so timesteps needed to be discarded. 

As initial data we choose a shell of matter described by a Gaussian scalar field: 



Its conjugate momentum is initially set to zero. This choice of initial data actually corre- 
sponds to a linear combination of an ingoing and outgoing pulse. As seen in Fig. ([T|, by the 
time the initial horizon forms the two pulses are well separated and the mass function and 
the lapse function are spatially constant to a close approximation. The spacetime between 
the pulses is therefore very nearly Schwarzschild. By fine tuning the initial data one can in 
principle consider a pulse that is entirely ingoing, but the presence of an outgoing pulse at 
a large distance from the black hole turns out to be useful to the present discussion. 

The parameters A, B and r are the initial parameters and determine the size of the 
outer apparent horizon. For the particular results we are presenting here we used the initial 
parameters (A, B, ro) = (0.185, 0.3, 3.0). This yields a black hole whose final mass (as 




(51) 




(52) 
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(a) Mass density (b) Lapse function 

FIG. 1: Plots of mass density and lapse function as a function of radius after the formation of 
the trapping horizon at time t = 1.73. The dotted lines are the two apparent horizons (ri n = 0.64 
and r ou t = 1.01). The lapse function has been normalized to unity at spatial infinity. The inner 
pulse is ingoing whereas the outer pulse is outgoing which is indicated with the arrows. They are 
sufficiently separated at this stage so that the region in between is approximately Schwarzschild. 
The mass function is constant in this region, while the lapse is spatially constant but increases 
with time. 



determined by the mass function between the pulses) is about 0.51. The ADM mass of the 
spacetime, including the outgoing pulse, is about 0.65. 

One of the main advantages of PG coordinates over null coordinates that are typically 
used for such simulations is that the slicings extend into the trapping region and the code 
can be run beyond horizon formation. Fig. (J2J) maps out the trapping horizon as a function 
of PG time. One can see that the horizon forms initially at about t=1.66, and evolves into 
an inner and outer apparent horizon. As mentioned above, the former necessarily moves 
inward while the outer horizon grows until it reaches its final value after all the ingoing 
matter has fallen through. This occurs by about t=1.78. It is not unusual for additional 
pairs of apparent horizons to form, as happens in the black hole interior in this evolution 
at t=1.73. Once the innermost horizon reaches the origin, the singularity forms and the 
code terminates. It can be verified that the trapping horizon for our simulation is spacelike 
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everywhere. 
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FIG. 2: Plot of location of outer apparent horizon (solid line) as a function of PG time. Note that 
between t = 1.73 and t = 1.79 there are four apparent horizons. The singularity forms and the 
code terminates when the innermost horizon hits the origin at t = 1.87. 



The lapse function has been normalized to unity at spatial infinity, i.e. outside the 
outermost pulse. As the system evolves, the lapse function in between the pulses can depend 
on both r and t. After the pulses are well separated, Birkhoff 's theorem requires the lapse 
between the pulses be spatially constant. However, the time dependence persists. Fig. 



(4(a) ) shows the lapse function on the outer apparent horizon as a function of PG time. 
It decreases rapidly while the horizon is growing. Once the matter has fallen through, it 
is spatially constant as expected but increases with time. This is of course a coordinate 
dependent, kinematical effect that is a consequence of the presence of an outgoing outer 
pulse. The effect can easily be illustrated analytically by considering the outgoing pulse to 
be infinitely thin and using the Israel junction conditions (see for example [31]). In this 
case, use of a single time coordinate on both sides of the shell while requiring continuity of 
the metric tangential to the null world line of the shell determines the discontinuity in the 
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lapse across the shell be given by: 



a in (t) = a out (t)(l - 2Af/ife(t)) (53) 

where -Rs(t) is the radius of the shell at time t. If a out = 1, then <j in asymptotes to one as 
the shell moves outwards. This is consistent with the observed increase in a on the outer 
horizon. In order to verify the details of the asymptotic behaviour one would have to run the 
code well past singularity formation. We have instead taken the simpler but equally effective 
route of doing a subcritical run in which all the matter reflects from origin and subsequently 
moves outward. Fig. ^ illustrates the time dependence of the lapse at a radius r b s near 
the origin. Although the lapse exhibits the expected asymptotic approach to unity once the 
shell has escaped (i.e. t > 3) the convergence is surprisingly slow. 

One can also choose a new time coordinate t in which the lapse is constant in the re- 



gion between the pulses. Fig. (4(b) ) plots the rescaled lapse ait) on the outer horizon in 
coordinates in terms of such a time coordinate. This corresponds to normalizing the local 
Killing vector d/dt to one in the region just outside the horizon but well inside the outgoing 



pulse. The fact that a deviates slightly from unity at large PG times in Fig. ( |4(b) ) is a 



consequence of the fact that not quite all of the matter has yet fallen through the horizon 
at singularity formation, so that a' is not precisely zero. 

The following graphs show the dynamical behaviour of the various definitions of surface 
gravity during the formation of the outer apparent horizon as it settles down to the final 
event horizon. 

Fig. ([5]) calculates the surface gravity using the rescaled lapse <r, i.e. from the point of 
view of an observer in the approximately Schwarzschild region between the two pulses who 
is presumably unaware of the existence of the outgoing pulse. From this perspective the 
spacetime settles down to that of an isolated Schwarzschild black hole. As expected all the 
definitions asymptotically converge to the same value of about 0.5, but the there is some 
variety in the rate of convergence. The definitions without time derivatives start at zero 
on formation and are practically indistinguishable because the lapse function at the horizon 
does not change much during the evolution, kpq on the other hand depends explicitly on 
the time rate of change of a, which is relatively large and negative during the dynamical 
phase. The resulting surface gravity is therefore negative. The remaining two definitions 
start off slightly higher than the Killing value, because M is positive, and settle down very 
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FIG. 3: Plot of the lapse function corresponding to subcritical collapse for initial parameters 
(A, B, ro) = (0.1, 0.3, 1.0) which yield an ADM mass of about 0.02. The lapse function is evaluated 
at the observation point of r f, s = 0.77. There are two pulses: an ingoing and an outgoing. The 
ingoing pulse falls below r b s , bounces and is reflected from the origin before escaping to spatial 
infinity just like the first outgoing pulse. The structure in the time dependence of a(t) between 
and about 2 PG time is due to the passing of the inner pulse. After about 3 PG time units, 
all matter is outside r b s so that the spatial derivative of the lapse is approximately zero and the 
spacetime is very nearly flat. 

rapidly. 

Fig. (j6l) on the other hand uses PG coordinates in which the lapse is normalized to unity 
at spatial infinity, i.e. outside the outgoing pulse, as demanded by the standard Killing 
definition. This does not change significantly the initial dynamical behaviour of the surface 
gravity, but it does affect the final equilibrium value. The Hayward and Nielsen- Visser 
definitions settle down to a different constant value because they are missing the overall 
factor of the lapse function which is not unity at the outer horizon even after all the matter 
has fallen through. As argued above one expects the lapse to asymptote at long times to 
one, so that the definitions will in this case ultimately agree, in distinction to what happens 
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(a) Lapse Normalized at Infinity (b) Lapse Normalized outside Horizon 

FIG. 4: A plot on the left shows the lapse function at the horizon as a function of PG time 
normalized to be unity at spatial infinity. The plot on the right shows the rescaled lapse function 
a ("New Lapse Function") at the horizon as a function of the rescaled time t ("New PG Time"). 



for static, dirty black holes. 



VI. SUMMARY AND CONCLUSION 



We have presented an analytic and numerical study of several definitions of dynamical 
surface gravity. Our results are summarized in Table ([!]). Only two of the definitions yield 
zero for extremal horizons. Moreover, as indicated in the last two columns, the different 
classes of definition exhibit qualitatively different behaviour in the dynamical region. The 
PG definition (which contains an explicit b in the definition) is negative on formation, 
which may be sufficient grounds to rule. The definitions with M start slightly above the 
equilibrium or final value, whereas the two definitions that had neither a nor M start at 
zero on formation and take longer to converge. 

As expected, all definitions agree for static horizons with vacuum exteriors. More gen- 
erally (i.e. non- vacuum exteriors), even in the static case one must choose between locality 
(independence of matter content outside the horizon) and agreement with the Killing def- 
inition for all black holes with a global timelike Killing vector. This was illustrated with 
a couple of specific examples of dirty black holes and more graphically in our numerical 
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FIG. 5: Plot of various definitions of dynamical surface gravity on the outer horizon vs. new PG 
time with lapse chosen so that the timelike Killing vector is normalized to unity in between the 
two shells of matter. This is essentially equivalent to ignoring the outgoing pulse and treating the 
ingoing shell as isolated. In this case the asymptotic values of all definitions agree. 

analysis. The non-local definitions behave differently depending on whether one normalizes 
the Killing vector in the approximately Schwarzschild region directly outside the horizon 
(i.e. between the shells) or at spatial infinity. One somewhat unexpected result is that 
when the Killing vector is normalized at infinity the non-local definitions produce values for 
the surface gravity that increase with time even after all the matter has fallen through the 
horizon. This effect was attributed to the presence of the outgoing shell of matter between 
the observer at infinity and the event horizon and further highlights the non-locality of the 
standard definition. 

One important question concerns the magnitude of the differences. Although we con- 
sidered only massless scalar field collapse, our simulations nonetheless provide hints about 
the relative magnitude of the differences between the various definitions in more realistic 
astrophysical situations. The massless scalar field dynamics is scale invariant in the sense 
that the amplitude of the initial pulse is dimensionless, while the units for the width B and 
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FIG. 6: Plot of various definitions of dynamical surface gravity on the outer horizon vs. PG time 
with lapse chosen so that the timelike Killing vector is normalized to unity at spatial infinity. In 
this case the non-locality of some of the definitions is apparent: the asymptotic values do not agree 
even after the horizon has reached its final value. 

location r can be specified arbitrarily. For the values of the parameters represented by 
Figs. 1 through 5, the mass of the shell is about 20% of the final mass of the black hole, 
independent of the chosen scale. Thus if one chooses the units of the initial parameters to 
be, say, ten kilometers, then the final horizon radius is about 5 km, corresponding to a solar 
mass black hole, while the mass of the outgoing shell is about 0.2 of a solar mass. The 
difference between the local and non-local values of the final surface gravity, as represented 
by the asymptotic lines in Fig. 5, is about 20% when the outgoing shell is about 4 horizon 
radii from the origin. Although the different values are expected to converge at very long 
PG times, the convergence is quite slow, as illustrated in Fig. 
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TABLE I: Comparison of main properties of various definitions of dynamical surface gravity con- 
sidered. The column labelled "value on formation" indicates whether the surface gravity starts off 
as zero, negative, or greater than the final (or equilibrium) value which is necessarily positive. In 
all cases studied the surface gravity is a monotonic function of PG time during the formation of 
the horizon. 

It is important to note that while several of the surface gravity definitions we have inves- 
tigated are motivated by and explicitly evaluated on apparent horizons, there is in principle 
nothing to stop us from evaluating the same quantities on null dynamical causal horizons 
such as the dynamical event horizon. In this case we would not expect any qualitative 
difference in our results. 

Our results are intended to provide guidance with respect to a useful choice of definition 
for the surface gravity of dynamical horizons. If one requires the surface gravity of extremal 
horizons to be zero, then there are only two viable candidates: the Hayward definition and 
the null definition. For the initial data we used, both these definitions exhibited qualitatively 
similar behaviour during the dynamical phase of the horizon evolution. One must also choose 
between locality of the definition, i.e. independence from the matter content outside the 
horizon, and universal agreement with the standard Killing definition for static black holes, 
including "dirty" black holes. In order to proceed one must know what physical information 
black hole surface gravity is meant to provide. For example, the Killing definition of surface 
gravity determines the tension in a massless string held by an observer at infinity suspending 
a weight just above the horizon. It is not surprising that such a measurement is sensitive to 
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the mass distribution between the horizon and the observer. 

In a more modern context, it is the connection with black hole thermodynamics that 
makes surface gravity important. Most derivations yield the temperature of the black hole 
as measured at infinity which is non-local in the sense given above. On the other hand, 
Hay ward et al. [10] have argued that in static cases at least, there exists a local definition 
of surface gravity that provides a local measure of the temperature of the horizon. 

Much work remains to be done before this important issue can be completely resolved. 
For example, in order to fully interpret dynamical surface gravity in terms of Hawking 
temperature, one needs a formalism for dealing quantitatively with the thermodynamics of 
dynamical horizons, i.e. the non-equilibrium thermodynamics of black holes. This is an 
even more fundamental and difficult problem. 
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